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2-color QCD (quantum chromodynamics with jV c = 2) at finite temperature T and chemical 
potential fi is revisited in the strong coupling limit on the lattice with staggered fermions. The 
phase structure in the space of T, fj,, and the quark mass m is elucidated with the use of the mean 
field approximation and the 1/d expansion (d being the number of spatial dimensions). We put 
special emphasis on the interplay among the chiral condensate (qq), the diquark condensate (qq), 
and the quark density (q q) in the T-\i-m space. Simple analytic formulae are also derived without 
assuming \x nor m being small. Qualitative comparisons are made between our results and those of 
' recent Monte-Carlo simulations in 2-color QCD. 

o : 
o ■ 

(NJ ; I. INTRODUCTION 

c : 

Physics of matter under high baryon density is one of the most challenging problems in quantum chromodynamics 
(QCD) both from technical and physical point of view. So far there have been proposed various novel phases which 
include the 3 P2 neutron superfluidity the pion and kaon Q condensations, the deconfmed quark matter Q, 

• the strange matter |5j, the color superconductivity |6|], the ferromagnetic quark ordering [3 and so on. 

Unfortunately, analyzing these phases from the first principle lattice simulations is retarded due to the complex 
fermion determinant for 3-color QCD at finite baryon chemical potential with quarks in the fundamental representa- 
tion. The situation is however different for 2-color QCD in which the fermion determinant can be made real and the 
\q ] Monte-Carlo simulations are attainable as was first demonstrated in Because of this reason, 2-color QCD provides 
an unique opportunity to compare various ideas at finite chemical potential with the results from lattice simulations. 
^sO I One of the major differences between 2-color QCD and 3-color QCD lies in the fact that the color-singlet baryons 
are bosons in the former. This implies that the ground state of the 2-color system at finite baryon density in the 
color-confined phase is an interacting boson system, i.e. a Bose liquid, although the quark Fermi liquid may be realized 
at high baryon density in the color-deconfincd phase. How this Bose liquid changes its character as a function of the 
,JZh 1 temperature T, the quark chemical potential /1, and the quark mass m is an interesting question by itself and may 
Oh, also give a hint to understand physics of the color superconductivity in 3-color QCD in which the crossover from the 
Bose-Einstein condensate of tightly bound quark pairs to the BCS type condensate of loosely bound Cooper pairs 
<D \ may take place @ . 

• In the present paper, we revisit the thermodynamics of the strong coupling limit of 2-color lattice QCD with 
staggered fermions with chiral and diquark condensates which was orig inally studied in [H GH 

Our main 

• , purpose is to analyze its phase structure and the interplay among the chiral condensate (qq) , the diquark condensate 
(qq), and the baryon density (q'q) as functions of T, /1, and m. This would give not only a useful guide to the actively 
pursued lattice QCD simulations of the same system in the weak coupling [la . Il4l Il5j but also give us physical insight 
into the Bose liquid together with the other knowledge from the instanton liquid model [16j , the random matrix model 
, the chiral perturbation theory 0, 0] , and the renormalization group |20| . 
This paper is organized as follows. In Sec. [HJ we introduce 2-color QCD with one-component staggered quarks 
(which correspond to quarks with 4 flavors in the continuum limit). By introducing the auxiliary fields a and A 
corresponding to the chiral and diquark condensates, we derive an effective free energy F e g[a, A] in the 1/d expansion 
(d being the number of spatial dimensions) and the mean field approximation. The resultant free energy is enough 
simple so that one can make analytic studies at least in the chiral limit m = with finite T and /1, and in the zero 
temperature limit T = with finite /1 and m. Sec. IIIII is devoted to such studies and useful relations of the critical 
temperature and chemical potential are derived. In Sec. II V I we make numerical analyses on the chiral and diquark 
condensates as well as the quark number density as functions of T, //, and m. The phase structure of 2-color QCD at 
strong coupling is clarified in the three dimensional (T, /x, m) space. In Sec. summary and concluding remarks are 
given. In Appendices El and [5] we give some technical details in deriving and analyzing the free energy. 



II. FORMULATION 



In this section, we derive an effective action for meson and diquark fields starting from the lattice action with the 
staggered fermion. The readers who are interested in the results in advance can skip the derivation and directly refer 
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to the expression of the free energy given in Eq. (|27|l with Eq. i|24|) . 
The action on the lattice is given by 



S[U, X ,x} = Sg[U}+S f [U,x,x] 



which consists of the gluonic part, 



S G [U] = ^J2 l 1 -^^ 11 ^)} > U lu ,(x) = U}(x)Ul(x + 1>)U v (x + ii)U lt (x), 
9 x,n,u ^ c J 



(1) 



(2) 



and the fermionic part with finite chemical potential |2l| . 



1 d 

Sf[U,x,x] = m^2x(x)x(x) + tj X^X/ 7 ^ {x( x )Uj(x)x(x + j) - x(x + (x)x(x)} 

X X j — 1 

+ \Y1 nix) {x(x)e»U (x) X (x + 6) - x(x + 0)U^x)c-^x(x)} 



(3) 



X stands for the quark field in the fundamental representation of the color SU(iV c ) group and t/ M is the SU(j?V c ) valued 
gauge link variable, g is the gauge coupling constant and d represents the number of spatial directions which takes 3 
in reality. Later we sometimes use a notation x = (t, x) in which r (x) represents the temporal (spatial) coordinate. 
rjo(x) and rjj(x) inherent in the staggered formalism are defined as 



m (x) = 1 



(x) = (-l)Ef =1 x i 



(4) 



/i in Eq. © is the quark chemical potential, while the temperature T — (aA^J^ 1 with a being the lattice spacing and 
N T being the number of temporal sites. We will write all the dimensionful quantities in unit of a and will not write 
a explicitly. Since we are interested in 2-color QCD in this paper, we take N c = 2 in the following. 

It is worth mentioning here on the symmetry breaking pattern of 2-color QCD formulated on the lattice. The 
staggered fermion with a single component, which corresponds to N{ = 4 in the continuum limit, has the global 
Uy(l) x Ua(1) symmetry for m = \i = [22]. Here Uy(l) (Ua(1)) corresponds to the conservation of the baryon 
number (axial charge). For N c = 2, Uy(l) x Ua(1) is graded to a larger symmetry U(2) [13j,|22j. This is because the 
color SU(2) group is pseudo-real and the action possesses Pauli-Giirsey's fermion— anti- fermion symmetry |23| . 

Introduction of a finite chemical potential /i explicitly breaks the U(2) symmetry down to the original 
Uy(l) x Ua(1). Further introduction of a finite quark mass m retains only the Uy(l) symmetry. Listed in Ta- 
bleware the symmetries realized in various circumstances and their breaking patterns. 





m = 


m 


n = 


U(2) 

broken to U(l) with 3 NG modes 


Uv(l) 
not broken 




U V (1) x U A (1) 
totally broken with 2 NG modes 


U V (1) 

totally broken with 1 NG mode 



TABLE I: Symmetry realized in the single component staggered-fermion action for iV c = 2. 
patterns and the number of Nambu-Goldstone (NG) modes are also listed. 



Possible symmetry breaking 



Before addressing the computational details to derive the effective action for the meson and diquark system from 
the original action Eq. ft}, we shall summarize the actual procedure which is similar to that of [111 ]: 

[Step 1] Strong coupling limit g — > oo is taken. Then the gluonic part of the action S'g]^] in Eq. J2J) vanishes, because 
it is inversely proportional to g 2 . Consequently the gauge field remains only in the fermionic part @. 

[Step 2] Large dimensional (1/d) expansion is employed in the spatial directions in order to facilitate the integration 
over the spatial link variable Uj. The temporal link variable Uq is left untouched at this Step and will be exactly 
integrated out later in Step 4. 

[Step 3] Bosonization is performed by introducing the auxiliary fields a for XX an d A for XX- Then the mean 
field approximation is adopted for the auxiliary fields. Namely, a and A are regarded as spatially uniform 
condensates. 



3 



[Step 4] Integration with respect to x, x, and Uo are accomplished exactly to result in an effective free energy written 
in terms of a (chiral condensate) and A (diquark condensate). 

Some comments are in order here to emphasize the importance of the exact integration over Uq without using the 
1/d expansion in Step 4. Suppose we performed the 1/d expansion not only in the spatial directions but also in the 
temporal direction to perform the Uj and Uo integrations. Then, the quarks would be totally replaced by the mesons 
and diquarks, and the temperature would not enter into the free energy within the mean field approximation. In other 
words, one needs to consider meson and diquark fluctuations to obtain thermal effects in such an approach. In our 
procedure, on the other hand, the thermal effect can be incorporated even in the mean field level by exempting the 
1/d expansion in making Uq integration. Similar observation has been made also in an exactly solvable many fermion 
system pEj . in which it was shown that the approach similar to ours has better convergence to the exact result. 

Let us now describe each procedure in some details with a particular emphasis on the importance of the interplay 
between the chiral and diquark condensates. 



A. Strong coupling limit and the 1/d expansion (Step 1 and 2) 

After taking the strong coupling limit g — * oo, the partition function is written as 

Z = J Vl Xl x}V{U }V[U 3 }e- s ^ u ^ . (5) 

Because x an d X are fermion fields, the Taylor expansion of e~ SF ^ U ' x '^ generates at most 2 2Nc terms on each site x. 
Thus the integration with respect to the link variable could be performed exactly in principle. Instead, we truncate 
such an expansion here up to the leading order in the 1 / d expansion. As explained soon below, the lowest order of the 
Taylor expansion gives the leading contribution in the 1/d expansion. Then we can integrate each term with respect 
to Uj at each site x. As we stated in Step 2 of the above summary, Uq is left untouched at this stage: 



d[J7j(x)] exp 



1 d 

-jj^Z^'fc) {x(x)U j (x)x(x + j) - x(x + j)U](x)x(x)} 



3=1 



d j d 

1 + g E X a {x) X a {x)x b (x + j) X h (x e abX a (x)x b (x)e cdX c (x + j)x d (x + j) + h. c. + 0(l/d) 



3=1 



3 = 1 



exp 



0(l/d), 



(6) 



M(x)V M (x, y)M{y) + £ B(x)V B (x, y)B(y) 
. v y 

where the Latin indices a, b, . . . are summed over in the color space. In integrating the link variable, we have utilized 
the formulae for the SU(2) group integration, 



J d[U] = 1 , J d[U] U ab ul d = ~6 ad 5 bc , J d[U] U ab U cd = ~e ac e bd . 



The mesonic composite M(x) and the baryonic composite B(x) are defined respectively as 
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M(x) = ±6 ab x a (x) X b (x), B(x) = K ab x a (x) X b (x), B(x) 
and their propagators in the spatial directions are given by 



e ab x a {x)x b (x) 



1 d 

V M (x,y) = V B {x,y) = t ^ [ S y,x 



(7) 



(8) 



(9) 



3 = 1 



The reason why the bilinear forms of M(x) and B(x) correspond to the leading order of the 1/d expansion can 
be understood by changing the normalization as Vu(x,y)/d — > Vm(x, y) and M[x) — ► M(x)/\/~d so that the mesonic 
propagator Vm(x,u) is 0(1) in the large d limit. As a result, the more M(x) is contained in the higher order terms 
in the Taylor expansion, the more 1/Vd is associated with it. The same argument holds exactly for the baryonic 
composite B(x) with N c — 2, where B(x) is composed of two quarks. For N c > 2, however, the baryonic contribution 
is of higher order in the 1/d expansion as compared with the mesonic one. For further details on the 1/d expansion 
for general N c , see also the discussion in 22]. 
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B. Bosonization and the mean field approximation (Step 3) 

The resultant action JfjJ) describes the nearest neighbor interaction between the mesonic composites and between the 
baryonic composites. Since they are the four-fermion interactions, one may linearize them with a standard Gaussian 
technique with the auxiliary fields a and A; 



exp 



Y,M(x)V M (x,y)M(y) 



x.y 



J V[a] exp 



-J2{a{x)V u (x,y)a(y) + 2a(x)V M (x,y)M{y)} 



x.y 



(10) 



and 



exp 



Y,B(x)V B (x,y)B(y) 



= I V[A] exp 



■^2{A*(x)V B (x,y)A(y) - A*(x)V B (x,y)B(y) - B(x)V B (x, y)A(y)} 



x.y 



From the above transformations, it is easy to show the relation, 

(a(x)) = - (M(x)) , (A(a:)) = (B(x)) , (A*(x)) = (B(x)) , 



(11) 



(12) 



where we have intentionally chosen the definition of the field a so that the sign of {cr(x)} becomes positive for m > 0. 

Now we take a mean field approximation. Namely, we replace the auxiliary fields cr(x) and A(x) by the constant 
condensates a and A and ignore any fluctuations around the condensates. It is obvious from Eq. Ijl2|l that a should 
be identified with the chiral condensate and A be the diquark condensate. Then we can write the partition function 

as 



Z= / V[U }V[ X ,x]e 



-S'[U ,x,x;<rA] 



(13) 



with 



'[U , x, XI ^ A ]=J2 m x(x)x(x) + ~ {x(x)e tt U (x)x(x + 6) - X {x + 0)U f o (x)e "x(*)} 



^ CT 2 + d| A j a +daM ^ _ ^ A * B(x) 



-AB(x) 



(14) 



Since S [Uq, x, x; c, A] is in a bilinear form with respect to the quark fields x an d X, we can integrate out them 
immediately. The integration with respect to Uo, which seems to be tough at first glance, turns out to be feasible as 
demonstrated in the next subsection. 



C. Integrations over Xt Xt and Uo (Step 4) 



In order to complete the remaining integrals, we adopt a particular gauge in which Uq(t, x) is diagonal and inde- 
pendent of t (often called the Polyakov gauge) , 

U (t,x) = diag { e i0 iW/ N -, d**W/"A , with 6i{£) = -6 2 {x) . (15) 

Also we make a partial Fourier transformation for the quark fields; 

1 Nt 1 Nt (m-i) 

X (r ) x) = -=^e ifc -x(m ) x), X (r,x) = -= ]T e - ik -^(m,x) , fc m = . (16) 

* T m— 1 * T m— 1 r 

The anti-periodic condition in the temporal direction, x( T + N T ,x) = — X (t,x), is satisfied owing to the fermionic 
Matsubara frequency k m . 
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Substituting Eqs. I|15[) and (|16(l into the action l|14fl and taking the summation over r, we reach the action in a 
Nambu-Gor'kov representation; 



(17) 



where 



X"(m,x) = [X a M\ 
\X a (m,x)l 



and 



G^{m,n;9(g)) 



-A* 5 m ,iV T -n+l e o6 



V 



M + i sin fc. 



9q(a) 



1/' 



M - i sin I fe 

)l 
2 



9a(x) 



lfl 



A S m ,N T -n+l^ab 



\ 



J 



(18) 



(19) 



The indices a and b run from 1 through 2 in the color space. M denotes the dynamical quark mass (to be distinguished 
from the mesonic composite field M(x)) defined by 



M = m + —a . 

We can perform the Grassmann integration over x and y as [24 



= ][[V Det KKn^Cx))] 



(20) 



(21) 



Det stands for the determinant with respect to the color indices and the Matsubara frequencies. The square root of 
the determinant may be simplified as 



v/DctfG- 1 ] = J] 



I ) |A| 2 + {M + isin (k m + 6/N T - iyu)} • {M - isin (k m + 9/N T + i/j,)} 



(22) 



where 9\ = —9i = 9 is substituted. The product with respect to m can be performed using a technique similar to that 
in the calculation of the free energy in finite-temperature field theory in the continuum. The details of the calculation 
is given in Appendix 1X1 The result turns out to be a rather simple form, 



'Det [G~£(m,n;6)] = (cos 9 + cosh N T E_) ■ (cos 9 + cosh N T E^ 
Here E± is the excitation energy of quasi-quarks, 



E ± =arccosh y (1 + M 2 ) cosh /1 + (d/2) 2 |A| 2 ± Msinh fi 



(23) 



(24) 



with the dynamical quark mass M defined in Eq. (|20|l . 

Finally, all we have to do is to integrate this resultant determinant with respect to Uq, or 9, to derive the effective 
action; 



S eff [^A] = -logZ = Wlr 2 + !|A' 2 



E lo s 



d9(x) . 



2 2 1 1 J ^ [J 2ir 

x - X 

where we have used the SU(2) Haar measure in the Polyakov gauge <|15[) . 



sin 2 9(x) ^Det [G~^(m,n;9(x))\ \ , (25) 



7 rr :: J{ /" li^W 



2n 



(26) 



e{x)=e 1 {x)=-e 2 {x) 
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The Uq integration projects out the color singlet pairing of quarks and anti-quarks among arbitrary excitations in the 
determinant. Thus what excites thermally is no longer single quarks, but color singlet mesons or diquarks. In this 
sense, the strong coupling limit in 2-color QCD inevitably leads to the boson system. 
Eqs. I)23|) . I|25|) . and l|26|) immediately yield the effective free energy; 

F eS [a, A] = Seff/(£*) = \<? + ^|A| 2 - Tlog{l + 4cosh(S + /T) • cosh(£_/T)} , (27) 

where we have rewritten N T in terms of the temperature T(— 1/N T ). 

Although the quasi-quark energy E± given in Eq. (|24|l takes a complicated form, E± is reduced to a simple expression 
in the naive continuum a — > with a being the lattice spacing. Assuming A/a, Aa, pa <C 1 in this limit, E± amounts 
to a familiar form in the continuous space-time, 

E ± ~ arccosh (Vl + M 2 +W+ (d/2) 2 |A| 2 ± Mp\ -» y/(M ± p) 2 + (d/2) 2 |A| 2 . (28) 

The quasi-quarks are static in this framework because of the mean field approximation. 



III. ANALYTIC RESULTS ON THE PHASE STRUCTURE 



A. Case in the Chiral Limit 



In the chiral limit m = with zero chemical potential p — 0, the free energy in the mean field approximation given 
in Eq. Ij27(l is a function only in terms of a 1 + | A| 2 . As a result, the free energy is invariant under the transformation 
mixing the chiral condensate with the diquark one. This corresponds to a subgroup of the U(2) symmetry of the 
original action at m — p = in Table [I] Because of this symmetry, the chiral condensate is indistinguishable from 
the diquark condensate for m = p = 0, so that a state with finite a can be arbitrarily transformed to a state with 
finite A. 

Finite p would act on the free energy as an external field tending to make the diquark condensation favored. Con- 
sequently even an infinitesimal introduction of p leads to the diquark condensation phase with zero chiral condensate 
in the chiral limit. This is a peculiar feature of 2-color QCD and is in a sharp contrast to the 3-color QCD. In fact, 
we can prove that the minimizing condition for our F e e[a, A] does not allow non-zero a at any T and p in the chiral 
limit. See Appendix iBl for the proof. 

Taking the fact that a vanishes in the chiral limit for granted, the free energy is simply written as 

d. 



F cff [A] = ^|A| 2 -Tlog{l + 4cosh 2 (£: /T)} , E = arccosh \J (d/2) 2 |A| 2 + cosh 2 fij . (29) 

Assuming that the finite T phase transition with fixed \x is of second order, which will be confirmed numerically later, 
the critical temperature T c may be determined by expanding _F c ff [A] in terms of |A| 2 and extracting the point where 
the coefficient of |A| 2 changes its sign. Since the expansion reads 



f„[a] . -Ti„ g{3 ^cosMWr,} + {I - 3 + 2 S^ m a 'tX ) } + (|A| 



one finds 



I / 3 sinh 2 2/i + d Ad 2 + 5 sinh 2 2u 

TJu) = 2u < arccosh ^ 

W \ \ 2d 2 -2 sinh 2 2/i 



(30) 



(31) 



The values of T c (/x) for two typical cases are 

6 1 

T c (u = 0) = -, TJfx = fj, c = - arcsinhd- 0.91) = (32) 
5 2 

for d = 3. T c is a monotonically decreasing function of /j, connecting the above two points. This will be discussed 
later in Sec. IIVI together with the case at finite m. 

It is interesting to calculate the quark number density p which is defined as 

^ _ 6>F cff _ 8 cosh (Eq/T) sinh (E /T) dE Q 
P ^ 1} - l + 4co S h 2 (E /T) ■ dp ' {66) 
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where Eq is the quasi-quark mass in the chiral limit defined in Eq. (|29|l . Although the above expression seems to be 
a little complicated, it can be simplified by using the gap equation; 

Q = aFeff = d A _ 8cosh(E /T)smh(E /T) _ (ri/2) 2 A _ dEo 
<9A* 2 1 + 4cosh 2 (T^o/T) ' 2cosh/^sinh^ ' dp' 

In general the gap equation may have two solutions, A = and A ^ 0. The former (the latter) is the solution to 
minimize the free energy for T > T c ([i) (T < T c (fi)). It can be alternatively said that the former (the latter) is the 
solution for \x > fi c (T) (/i < /i c (T)) with the critical chemical potential /i c as a solution of Eq. (|31|1 in terms of /x. 

By eliminating most of the complicated part of Eq. (|33|1 by means of the gap equation (PU, we have the following 
expressions; 

{(2 sinh 2/i) /d for fj, < /z (T) (A ^ 0) , 
4sinh(2 / VT) (35) 
3 + 2co S h(2 M /T) f ° r (A = 0) - 

At zero temperature, Eq. (|33|) is reduced to 

|(2sinh2 M )/d for M < Mc (T = 0) , 
^ T = 0) = { 2 for M >Mc(T = 0) (36) 

with 

iz c (T = 0) = iarcsinhd, (37) 

at which the quark number density p gets to be saturated to two. Note that quarks on each lattice site have only two 
degrees of freedom in the color space because the spin and flavor degrees of freedom are dispersed on the staggered 
lattice sites. Therefore "two" is the maximum number of quarks placed on each lattice site due to the Pauli exclusion 
principle. 



B. Case at Zero Temperature with Finite m 

Even a small quark mass could modify the phase structure substantially from that given in Sec. IIII Al This is 
because the pion, that is the Nambu-Goldstone mode associated with the chiral Ua(1) symmetry breaking, comes to 
have a finite mass, oc m 1 ^ 2 . As a result, as far as the quark chemical potential is smaller than a threshold value, 
the vacuum is empty and the diquark condensate vanishes; p(/i, T = 0) = and A.(jx, T = 0) = 0. 

To clarify such threshold effects in more detail, let us now focus our attention on the T = system with finite m. 
We derive analytic formulae for lower-critical chemical potential fi l ° w at which A starts to be non-vanishing and the 
upper-critical chemical potential /i" p at which A ceases to be non-vanishing. 

The free energy at T = with finite m takes a simple form; 

F eS [a,A] = A -a 2 + ^|A| 2 - (E + + E_) . (38) 

The expansion of F c ff[<r, A] in terms of |A| 2 near the threshold gives a condition to determine the critical chemical 
potential /i c , 

cl d 2 ( 1 1 

2 8 cosh E cosh fi c \smh(E + fi c ) \sinh(E — fi c )\ 

where E = arccosh (VT+ M 2 ) and M — m + da/2 as defined in Eq. (|20|l . In Eq. (|39|l the dynamical mass M, or a, 
is determined by the condition to minimize the free energy at the threshold. We can reduce that free energy much 
more by putting A = 0, 




FM = ^a 2 -(E + ii+\E-p\) 



(40) 
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The stationary condition of the free energy F c g [a] with respect to a yields the following chiral gap equations; 

2 {{l + M 2 )- 1 ' 2 for n<E, 

-JM-m) = < r (41) 

cT ' | for p>E. y ' 

Combining Eq. JSHJ and Eq. ijlTjl , we find 



/4° w = arccoshVl + mM for p < E , 



Mc ' -i^ p = arccosh\/l + if 2 for /x > J5 , ^ 



where we define K as the solution of the equation, 

^-^) = (i + ^ 2 r 1/2 . (43) 

With finite m at T = 0, the empty vacuum gives p = and A = as long as p < p l ° w . On the other hand, the 
non- vanishing value of A is possible for p]? w < p < /i" p . For p > p^ p , the saturation of the quark number density 
occurs leading to p = 2 and A = 0. These behaviors will be confirmed numerically in Sec. II VI 

Note that we have not assumed the quark mass m to be small, which is in contrast to the approaches based on 
chiral perturbation theory [Isl IT^|. Therefore, Eq. 1421) can relate the critical chemical potentials to arbitrary values 
of m. For sufficiently small m, we have a relation; 

p}° w = m 1 / 2 ■ l { ' 2 ~\ , Mc P = ^arcsinhd. (44) 

The former may be rewritten as p}° w — m-^jl with obtained from the excitation spectrum in the vacuum up 
to the leading order of the l/d expansion. This observation is consistent with the discussion given in 0|. The latter 
relation is simply rewritten as p,^ p = p. c {T = 0) with p c {T = 0) defined in Ea. (|37Jl . 
For sufficiently large m, we find asymptotically, 

/^ ow ~ arccoshm. (45) 

This is because M ~ K ~ m for large m as is evident from Eqs. 141(1 and l|4H|l. 



IV. NUMERICAL RESULTS ON THE PHASE STRUCTURE 



In this section, we determine the chiral condensate <r and the diquark condensate A numerically by minimizing the 
effective free energy in Eq. 1)2 7[l. The quark number density p is also calculated numerically. The results are shown 
in Figs, ^and [21 for m — and m = 0.02, respectively. The phase diagrams of the system in the T-fi plane, in the 
p-m plane, and in the three dimensional T-p-m space are also shown in Figs. [31 El and [3 respectively. 



A. Case for m = 



Let us first consider the diquark condensate and the quark number density as functions of p and T shown in the four 
panels of Fig.^ (Note that the chiral condensate is always zero in the chiral limit as we have proved in Appendix iBl) 

At T — (the upper left panel), the diquark condensate A decreases monotonously as a function of p and shows 
a second order transition when p becomes of order unity. On the other hand, the quark number density p increases 
linearly for small p and grows more rapidly for large p until the saturation point where quarks occupy the maximally 
allowed configurations by the Fermi statistics. (See the analytic formula given in Eq. 1(36(1 .) Those behaviors of A 
and p are also observed in the recent Monte-Carlo simulations of 2-color QCD [lj|. Note that the rapid increase of p 
near the upper-critical chemical potential p^ takes place even in the strong coupling limit as shown here. Namely it 
does not necessarily be an indication of the existence of free quarks at high density unlike the suggestion given in the 
last reference in |l4| . 

As T increases, the magnitude of the diquark condensate decreases by the thermal excitations of quark and anti- 
quark pairs in the last term of Eq. ((27(1 . which is shown in the lower left panel of Fig.^ It is worth mentioning here 
that the diquark condensate disappears even before the complete saturation (p = 2) takes place. 

Next we consider the diquark condensate as a functions of T for two typical values of the chemical potential in the 
right panels of Fig.^ The diquark condensate shows a second order transition at T c given analytically by Eq. ((31|) . 



9 




FIG. 1: Chiral condensate a, diquark condensate A, and the quark number density p for m = with d = 3. In the left panels 
they are plotted as functions of the chemical potential p for two typical values of temperature. In the right panels, they are 
plotted as functions of temperature T for two typical values of p. All the dimensionful quantities are in unit of the lattice 
spacing a, which is implicitly understood in all the figures in this section. 



B. Case for ra/0 

In the upper left panel of Fig. [2J the chiral condensate, the diquark condensate, and the quark number density are 
shown as functions of p for small quark mass m — 0.02 at T = 0. As we have discussed in Sec. IIII Bl there exists a 
lower-critical chemical potential p}° w given by Eq. (|42|) . Both A and p start to take finite values only for p, > p}° w at 
T = 0. 

One can view the behavior of chiral and diquark condensates with finite quark mass as the manifestation of two 
different mechanisms: One is a continuous "rotation" from the chiral direction to the diquark direction near p = p}° w 

or T — T c (p) with \/ a 2 + A 2 (shown by the dotted line) varying smoothly. The other is the saturation effect which 
forces the diquark condensate to decrease and disappear for large p, as seen in the previous case of m = 0. 

The "rotation" can be understood as follows: The free energy in the mean field approximation at small m and 
p, has an approximate symmetry which mixes the chiral condensate with the diquark one as we have discussed in 
Sec. IIII A1 The effect of m (p) is to break this symmetry in the direction of the chiral (diquark) condensation favored. 
Therefore, a relatively large chiral condensate predominantly appears for small p region. 1 Just above the lower- 
critical chemical potential p}° m , the chiral condensate decreases, while the diquark condensate increases because the 
effect of p, surpasses that of m. As p becomes large, the diquark condensate begins to decrease in turn by the effect 
of the saturation and eventually disappears when p exceeds the upper-critical value p^ p (order of unity for T = 0) . 

Similar "rotation" and the saturation effect are also seen at finite T as shown in the lower left panel of Fig. At 
finite T, both the chiral and diquark condensates are suppressed due to the effect of temperature and the diquark 
condensate disappears before the complete saturation occurs. 



1 Note that the chiral symmetry is explicitly broken by m, thus the chiral condensate is always non-vanishing although it is suppressed 
in magnitude at high T and /i. 
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FIG. 2: Chiral condensate a, diquark condensate A, and quark number density p for m — 0.02 with d = 3. All the definitions 
are the same with those of Fig. except that the dotted line indicates a total magnitude of the condensates \J a 2 + A 2 . 



Next we consider the chiral and diquark condensates as functions of T (the right panels of Fig. 0). At low T, 
both the chiral and diquark condensates have finite values for p = 0.2 and 0.5. The diquark condensate decreases 
monotonously as T increases and shows a second order transition. On the other hand, the chiral condensate increases 
as the diquark condensate decreases so that the total condensate \/ o 2 + A 2 is a smoothly varying function of T. 
The understanding based on the chiral-diquark mixing symmetry is thus valid. An interesting observation is that the 
chiral condensate, although it is a continuous function of T, has a cusp shape associated with the phase transition of 
the diquark condensate. 

Finally let us compare the to = case in Fig. ^ and the to = 0.02 case in Fig. [21 Looking into two figures 
at the same temperature or chemical potential, we find that the diquark condensate A for to = and the total 
condensate \/ a 1 + A 2 for to = 0.02 have almost the same behavior. This indicates that although the current quark 
mass suppresses the diquark condensate, the price to pay is to increase the chiral condensate so as to make the total 
condensate insensitive to the presence of small quark mass. Restating this by use of the radial and the angle variables 
defined by Asin0 = a and Acos0 = |A| as in Eq. HB1|) . a small to hardly changes the behavior of A but shifts </> from 
zero. 



C. Phase Diagrams 



Now we show the phase diagram of the strong coupling 2-color QCD in the T-p plane in Fig. [21 The solid line 
denotes a critical line separating diquark superfluid phase A / and the normal phase A = in the chiral limit 
to = 0, which is determined analytically by Eq. (|31[) . The dashed line denotes the critical line determined numerically 
for to = 0.02. The phase transition is of second order on these critical lines. The chiral condensate a is everywhere 
zero for to = 0, while it is everywhere finite for to ^ 0. In the latter case, however, a is particularly large in the region 
between the solid line and the dashed line. Shown in Fig. 0] is a phase diagram in /i-to plane at T = 0. The lower 
right of the figure corresponds to the vacuum with no baryon number present, p = 0. On the other hand, the upper 
left of the figure corresponds to the saturated system, p = 2, in which every lattice site is occupied by two quarks. 
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FIG. 3: Phase diagram of the strong coupling 2-color QCD 
in the T-/i plane. Solid (dashed) line denotes the critical line 
for diquark condensation for m = (m = 0.02). 
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FIG. 4: Phase diagram of strong coupling 2-color QCD in 
the p,-m plane at T = 0. Two solid lines separate the region 
where A = from the region where A / 0. 



There is a region with < p < 2 and A ^ bounded by the above two limiting cases, which is separated by two 
critical lines given in Ea.i|42p. 

Finally, bringing all the discussions together, the phase structure in the three dimensional T-fi-m space is shown in 
Fig. [SJ The diquark condensate has a none- vanishing value inside the critical surface and the phase transition is of 
second order everywhere on this critical surface. The second order phase transition is consistent with other analyses 
employing the mean field approximation; the random matrix model and the Nambu— Jona-Lasinio model |17| . and 
also the chiral perturbation theory at T = 18]. On the other hand, Monte-Carlo simulations of 2-color QCD [T^j 
show indication of a tricritical point in the T-fi plane at which the property of the critical line changes from the second 
order to the first order as /i increases. This is also supported by the chiral perturbation theory beyond the mean field 
approximation at finite T [l9| . 

Aside from the fact that we are working in the strong coupling limit (g — ► oo), though the lattice simulations are 
aiming at the weak coupling limit (g — *• 0), it would be of great interests to go beyond the mean field approximation 
in our analysis and study corrections to the phase structure given in Fig. [SJ 




FIG. 5: Phase structure of strong coupling 2-color QCD in the T-^i-m space. The surface represents the critical surface for the 
diquark condensation, which separates the region where A = from the region where A ^ 0. 
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V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have studied the phase structure of 2-color QCD in the strong coupling limit formulated on the 
lattice. We have employed the 1/d expansion (only in the spatial direction) and the mean field approximation to 
derive the free energy written in terms of the chiral condensate a and the diquark condensate A at finite T, p, and m. 
A major advantage of our approach is that we can derive rather simple formulae in an analytic way for the critical 
temperature (T c ), the critical chemical potentials (p l ° w and /i" p ), and the quark number density (p), which are useful 
to have physical insight into the problem. Although our results are in principle limited in the strong coupling, the 
behavior of cr, A, and p in the T-p-m space has remarkable qualitative agreement with the recent lattice data. Since 
we do not have to rely on any assumption of small p nor small m in our approach, the strong coupling analysis 
presented in this paper is complementary to the results of the chiral perturbation theory and directly provides a 
useful guide to the lattice 2color-QCD simulations. 

There are several directions worth to be explored starting from the present work. Among others, the derivation 
of an effective action written not only with the chiral and diquark fields but also with the Polyakov loop, L(x), is 
interesting because this will allow us a full comparison for all the physical quantities calculated in our formalism and 
lattice simulations [l4|. This is a straightforward generalization of the works in j2^, which allows us to analyze (L(x)) 
under the influence of chiral and diquark condensates. 

Also calculation of the meson and diquark spectra including pseudo-scalar and pseudo-diquark channels will give us 
deeper understanding about the non-static feature of the present system. Another direction is to study the response 
of the Bose liquid discussed in this paper under external fields. Although mesons and diquarks are color neutral in 
2-color QCD, they can have electric charge so that the external electromagnetic field leads to a non-trivial change 
of the phase structure if the field intensity is enough strong. Generalization of the whole machinery in this paper to 
strong coupling 3-color QCD is also one of the challenging problems to be studied in the future. 
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APPENDIX A: SUMMATION OVER THE MATSUBARA FREQUENCY 



Let us first take the logarithm of Eq. (|2"2l ; 



log y/BetiG- 1 ] = lo § [- cos2 ( fc ™ + °/ N r) + 2Y C0S (*m + 9/N T ) + X 2 ] 

m— 1 

Here X and Y arc defined as 



X 2 = cosh 2 \i + M 2 



Y = M sinh /i . 



(Al) 



(A2) 



By differentiating log y/Det[G 1 ] with respect to X, we obtain 







X 



v rn— 1 



E 



i 



cos (k m + 9/N T ) -Y + \/X 2 + Y 2 cos (k m + 9/N T ) - Y - \/X 2 + Y 2 



(A3) 



Because cos(fc m + 8/N T ) is invariant under the shift to — ► m + N T , we can make the summation of m over the range 
to = — oo to to — +oo with an appropriate degeneracy factor £1. Then the residue theorem enables us to replace the 
summation by a complex integral as 



d 

dX 



log y/Dct[G- 



X 1 
VX 2 + Y 2 ~& 



1 



27ri cos(z + 9/N T ) + fc_ 



-iN T x - 1 ~iN T 
1 + e iN - z ^ - sin z 1 + e iN ^ 



dw 

2ni cos(w + 6/Nt) - k+ 



1 -iN T ^ 1 -iiV T 

l/N T )-k+ l + e iN - w + 2s - s [nw 1 + e iN - 1 



(A4) 



where 



fc± = VA^ + y2 ± y 



(A5) 



Owing to the infinite range of the summation over to, we can choose the closed contour at infinity for the complex 
integrals with respect to z and w, and thus such complex integrals go to zero, z and w are the residues satisfying 



cos (z + 9/N T ) + = , cos (w + 6/N T ) -k+=0. 
Solving these equations, we obtain 

Z + 9/N T = ±lE- +7T + 27T77, W + 8 / N T = ±\E + + 27T72, (?7, = — OO, ■ ■ • , Oo) . 

where 

E± = arccoshfc± . 

Substituting Eq. HA7|> into Eq. (|A4|) . we obtain 
d 



(A6) 
(A7) 
(A8) 



dX 



log A/Det[G- 



X 



VX 2 + Y 2 V 



OO 

5 E 



71 — — OO 



1 



-U k 2 _ -1 



-iiV T 



-iAL 



]^ _|_ c -A r T £;_-i6(+iA r T 7r } _|_ C N T E- -i8+iN T -K 



1 



i^fci - 1 



-iiV T 



-iiV r 



9£ 4 



M \l + c - N - E -- ie l+e N - E -- ie J dX \l + e - N - E +- w l+e^^- 1 
A [bg { e ^ B - + e~ ie } + log {e~ N ^ E - + e- w } + log {e^+ + e- ie } + log {e"^ E + + e - i9 }] 

— [- 2i6» + 2 log 2 + log {cos 9 + coshN T E_} + log {cos 9 + coshN T E + }} . 



(A9) 
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We note that the degeneracy factor fl is just canceled by the infinite degeneracy of the summation on n. Also we 
have used the fact that N T must be an even integer for the staggered fermion. After the integration with respect to 
X , we find that -\/Det[G -1 ] is expressed in a rather simple form up to irrelevant factors, 

^Det [G~^(m,n;9)] = (cos<9 + cosh N T E_) ■ (cos(9 + cosh N T E+) . (A10) 
This is the final form shown in Eq. (12MI) in the text. 



APPENDIX B: PROOF OF a = AT THE GLOBAL MINIMUM OF F eB [a,A]\ m=0 

For to = 0, it is useful to rewrite the free energy Eq. (|27|) in terms of the radial and angle variables, 

cr = Asin0, |A|=Acos<£. (Bl) 

Then the free energy is 

F cS [\,4>] = ^A 2 -Tlog{l + 4cosh(£ + /T)-cosh(£_/T)} , (B2) 
with the quasi-quark energy given by 



E± = arccosh (d/2) 2 A 2 +cosh 2 /j+ (rf/2) 2 A 2 sinhVsin 2 0± -Asinh/jsin^J . (B3) 

What we are going to prove here is that F c g [A, 4>\ has the global minimum at <j) = for arbitrary T, /i, and A. Actually 
we can prove the following statement in a more abstract expression; 

cosh{aarccosh(&e c )} • cosh{aarccosh(6e _c )} has the global maximum at c = , 
where a > 1 , b > 1 , and — log b < c < log b . 

We can apply this corollary to our problem to prove that <f> = is the global minimum of the free energy. Obviously 
4> = turns out to be the global maximum of cosh(_E + /T) ■ cosh(£'_/T) once we substitute, 



a = l/T, b = y / (d/2) 2 A 2 + cosh 2 /j, c = log \Jl + (d\/2b) 2 sinh 2 /xsin 2 <j> + (dX/2b) sinh^sin0 



(B4) 



Thus we have proven that F c h[\, <f>] is globally minimized at <fi = 0, in other words, the chiral condensate vanishes in 
the chiral limit. 
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